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Neutrino flavor oscillations in the presence of ambient neutrinos is nonlinear in nature which leads 
to interesting phenomenology that has not been well understood. It was recently shown that, in 
the two-dimensional, two-beam neutrino Line model, the inhomogeneous neutrino oscillation modes 
on small distance scales can become unstable at larger neutrino densities than the homogeneous 
mode does. We develop a numerical code to solve neutrino oscillations in the multi-angle/beam 
Line model with a continuous neutrino angular distribution. We show that the inhomogeneous 
oscillation modes can occur at even higher neutrino densities in the multi-angle model than in the 
two-beam model. We also find that the inhomogeneous modes on sufficiently small scales can be 
unstable at smaller neutrino densities with ambient matter than without, although a larger matter 
density does shift the instability region of the homogeneous mode to higher neutrino densities in the 
Line model as it does in the one-dimensional supernova Bulb model. Our results suggest that the 
inhomogeneous neutrino oscillation modes can be difficult to treat numerically because the problem 
of spurious oscillations becomes more severe for oscillations on smaller scales. 

PACS numbers: 14.60.Pq 


I. INTRODUCTION 

The observation of neutrino flavor oscillations has established that neutrinos have non-vanishing masses and that 
their propagation (or mass) eigenstates are linear combinations of weak-interaction states. One of the remarkable 
successes in the field of experimental particle physics has been the measurement of all neutrino mixing parameters 
except the sign of the atmospheric mass-squared difference and the CV violation phase. 

Neutrino flavor oscillations can be modified by the potential due to the presence of electrons and nucleons in the 
medium or the presence of ambient neutrinos 00. There are two major differences in the phenomenology of 

neutrino flavor oscillations due to the presence of ambient neutrinos as opposed to ordinary matter. Firstly, unlike 
ordinary matter the presence of ambient neutrinos makes the neutrino flavor evolution nonlinear in nature. Secondly, 
in the case of ordinary matter the electrons and nucleons are usually non-relativistic, and the potential experienced 
by neutrinos is independent of direction to a very good approximation (see [6] for a review). This is not true in the 
case of neutrino-neutrino self-interaction. These two differences make neutrino oscillations in the neutrino medium 
very interesting and at the same time a challenging problem solve. It has been shown that a dense neutrino gas can 
undergo flavor oscillations collectively 00 . 

The effects of neutrino-neutrino interaction can be important in extreme environments with large neutrino densities 
like that in the interior of a core-collapse supernova. It was discovered in numerical simulations that the neutrino flavor 
evolution can be dramatically different for the normal and inverted hierarchies inside supernovae mm- However, 
in order to make the numerical simulations manageable, a simplified one-dimensional supernova model called the 
(neutrino) Bulb model was used. There are several effects that are not taken into account in the Bulb model, although 
they can modify neutrino oscillations significantly. For example, it has been found that the neutrino emission with 
an (approximate) axial symmetry around the radial direction can evolve into a configuration with a large axial 
asymmetry HUE]. It has also been suggested that the back-scattering of neutrinos from the nucleons in the envelope 
of the supernova can lead to significant modification of the neutrino potential m- 

In this paper we investigate the physics of collective neutrino oscillations in the neutrino Line model with two 
spatial dimensions. The study of this model can provide us with useful insights into the qualitative differences of the 
phenomenology of collective neutrino oscillations in models with one and multiple spatial dimensions. This study is a 
generalization of the work done for the two-beam Line model in Ref. M where only two neutrino beams are emitted 
from each neutrino source point. 
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FIG. 1. A schematic diagram of the two dimensional (neutrino) Line model. Each point on the a>axis or the “neutrino Line” 
emits neutrino beams with emission angles d within range [—'dmax, i?max]- 

II. THE NEUTRINO LINE MODEL 
A. Equations of motion 

In the stationary, two-dimensional (neutrino) Line model neutrinos and antineutrinos are emitted from the x-axis 
or the “neutrino Line” and propagate in the x-z plane (see Fig. [lj. We assume that the neutrinos and antineutrinos 
are of single energy E and the same normalized angular distribution g(d) such that the number fluxes of the neutrino 
and antineutrino within angle range [$, d + dd] are riyg{d)dd and nyg(d)dd, respectively, where d is the emission 
angle of the neutrino beam, and n v and rip are the (constant) total number densities of the neutrino and antineutrino, 
respectively. The flavor quantum states of the neutrino and antineutrino of emission angle d and at position (x, z) 
are given by density matrices p$(x,z) and p$(x,z), respectively [15]. We use the normalization condition 

tip = tip - 1 (1) 

such that the diagonal elements of a density matrix give the probabilities for the neutrino or antineutrino to be in the 
corresponding weak-interaction states. With these conventions the self-interaction potential for p$(x,z) in the Line 
model can be written as 


H 





[pw (x, z) - apw (x, z)\ [1 - cos(d - d')]g(d') dd f , 


( 2 ) 


where p = V2Gf^u with Gf being the Fermi coupling constant, and a = n„ jn v . In the Line model the strength of 
the neutrino self-interaction p is constant. In realistic astrophysical environments such as core-collapse supernovae, 
however, p can decrease with increasing distance from the neutrino source. 

The flavor evolution of the neutrino or antineutrino is governed by the equation of motion (EoM) 

• V p $ — [H vac T H mat + H (3) 


where is the unit vector that denotes the propagation direction of the neutrino with emission angle $, and H vac 
and H mat are the standard vacuum mixing Hamiltonian and matter potential, respectively. In this work we assume 
the mixing between two active neutrino flavors v e and v T with small vacuum mixing angle # v 1. Therefore, 


= H, 


^mat 


(A - rjLu) (\ 0 
2 VO -1 


(A - tjuj) 


-<? 3 , 


( 4 ) 


where A = \/2G Fn f . with n e being the net electron number density, rj is a parameter which takes a value of either +1 
or —1 for the normal neutrino mass hierarchy (NH, the mass-squared difference Am 2 > 0) or the inverted hierarchy 
(IH, Am 2 < 0), and uj = |Am 2 |/2£’ is the vacuum oscillation frequency of the neutrino. Eq. § can also be written 
in a more explicit form: 


i^osddz + sin dd x )p$ 




- cos(i? - #')] [P0> 


OCP0>, pv\ d-d'. 


( 5 ) 


The EoM for p$ is the same as Eq. <§ except with replacement uj —c a. 
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As in Ref. na we impose a periodic boundary condition along the x axis such that p$(x + L,z) = p$(x,z) and 
p$(x + L,z) = pft(x, z). It is convenient to recast the x-dependence of the neutrino density matrix in terms of Fourier 
moments: 


1 

Pm,o(z) = ~jr J e~ lkmX pi){x,z)dx, 


1 

Pm,i){z) = -J e~ lkmX p#{x,z)dx, 


where k m = 2tt m/L. It is straightforward to derive the EoM in the moment basis which are 

(A - 7]UJ) 


i cos Pd z p m ^ = fc m sin#p m> 0 + 




i COS 'ddzpmrf — km sin Ppm,ft H - 


+ p ^ ^ J [Pm',ft' &Pm',ft'iPm—m',ft\ [1 COs(P P )]p(P ) dP , 

m' 

(A + rjco) 


' [^" 3 > Pin , i )] 


/ \P” 

m f 


[1 - cos(i? - •&')]g(d ‘) d#'. 


(6) 


(7a) 


(7b) 


B. Collective modes in the linear regime 


We assume that the neutrinos and antineutrinos are emitted from the Line source in the electron flavor only. In 
the regime where neutrino oscillations are insignificant, the neutrino density matrices have the form 


P#(x,z) 


1 e# 

& 0 


p#(x,z) 


1 £0 
Tk 0 


(8) 


When there is a flavor instability, the off-diagonal elements of the density matrices grow exponentially, which can result 
in collective neutrino oscillations. In this section we apply the method of flavor stability analysis to the multi-angle 
Line model which was first developed in Ref. m- 
In the moment basis we have 


Pm, ft (^) 


^0 ,m 

r* 

— m,ft 


^m,ft 

0 


Pm,ft (^) 


do,m 

F* 

~—m,ft 


£m,ft 

0 


Keeping only the terms up to 0(e) in Eq. 0 we obtain 

icosPd z € m ^ = [A; m sin# + A - cjt? + (1 - a)/x^]e m>1? - p J [1 - cos(i? - 'd / )K e m,tf' - ae rn ^)g(P')dP', 
i COS Pd z €rn,ft — \J^m sin P H - A UJp ~\~ (1 C%) Pft\tm,ft P j [1 COS (P P )] (g rn,ft' m,ft ' ^)p{P ^)dP , 


where 


Pft = p J [l — cos($ — P')]g(P') dP' 


( 9 ) 


( 10 a) 

( 10 b) 


( 11 ) 


is the effective strength of neutrino self-interaction for the neutrino beam with emission angle P. As in the two-beam 
model, the flavor evolution of the neutrino fluxes in different moments is decoupled in the linear regime, although the 
evolution of the neutrino moments with different emission angles P are still coupled. 

Assuming that the mth neutrino moment oscillates with collective oscillation frequency we can write 


^m,i?(^) Qm,ft& m •> €m,ft(z') Qm,ft& m •> 

where Qm,ft and Qm,ft are ^-independent. Applying this ansatz to Eq. (JI(j|) we obtain 


( 12 ) 


IPm(pJ , P^)Qm,ft — ( p>m C m COS $ 8 m sin'd)p, 

Dm ( , P)Qm,ft — (^m C m COS P S rn sin P^)p 


(13a) 

(13b) 
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or 


where 


Q 


m 


Q 


m 


(g m - c m cos - s m sin d) g, 
D m (u,D) 

(dm - Cm cos - s TO sin #)/x 

D m (-w,i?) 


and 


Z} m (d=u;, d) = — Ucos'd + k m sind + A =p ujt] + (1 — <a)/i#, 


n m 

^m 


^(Qm^' )d(^ )dd , 

J (Qm,t?' -aQm,t)') c ost?'s(^)^ , . 
J(Qrn,v> - aQ m ^>)sind'g(d')dd'. 


Substituting Eq. (14) in Eq. 


( p~6] ) we obtain a characteristic equation for (a m ,c m ,s m ): 


( ^m[l] 1 Im [cos'd] / m [sind] \ / n m \ 

7 m [cosd] —7 m [cos 2 $] — 1 -7 m [cost? sin d] c m = 0, 

7 m [sin i?] —I m [cos d sin #] -7 m [sin 2 i?] - 1/ \s m J 


where 


imifm = I /WflW 

for arbitrary funciton /(d). Eq. © holds only when 


a/i 


_D m (u,d) D m (-u>,d)_ 


dd 


det 


7m [1] - 1 

I m [COS 1?] 

7 m [sin i?] 


—7 m [cos d] 

—7 m [cos 2 id] - 1 
— I m [cos 0 sin d\ 


-7 m [sin id] 
—7 m [cos'd sin d\ 
—Im [sin 2 i?] - 1 


= 0. 


(14a) 

(14b) 

(15) 

(16a) 

(16b) 

(16c) 

(17) 

(18) 

(19) 


For given to, A and g one can find a set of Qm (A, //) (i = 1,2,.,.) which satisfy Eq. (19) and which are the frequencies 
of the corresponding normal modes of collective neutrino oscillations. When 


K m = Im <Vm) 


( 20 ) 


is positive, the corresponding normal mode is unstable and its amplitude grows exponentially. If there exist multiple 
unstable modes, the mode with the largest exponential growth rate, 

C x =ma X (/ t W), (21) 


will eventually dominate. 


III. RESULTS IN THE LINEAR REGIME 
A. Numerical computation 

We develop a computer code to solve Eq. 0 numerically. In this code the continuous range ofd is represented as 
N discrete angle bins with central value d^ (i = 1,..., N ) and equal interval Ad. For an arbitrary function /(d) one 
has 


N 

/ /(d)dd —> Ad E m)- 

i =1 


(22) 
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FIG. 2. The evolution of |e m ,#|, the amplitudes of the off-diagonal elements of the neutrino moment matrices p m ,#(;z), in terms 
of propagation distance z for the inverted (left) and normal (right) neutrino mass hierarchies. The thick curves represent the 
numerical solution to Eq. 0 with Oth and 1000th moments only. The thin solid lines represent the exponential growth functions 
~ exp(ft™ ax z) predicted by the linear stability analysis. In these calculations we used the parameters listed in Eq. ( |24| , and we 
took the matter potential A = 0 and neutrino potential /jl/uj = 1500 (left) and 3000 (right) which is measured in the vacuum 
neutrino oscillation frequency uj. 


In this work we focus on the neutrino oscillations in the linear regime and the cases with a simple angular distribution 
which has isotropic neutrino fluxes within range [-^max,^max], i.e. 


aW = 


f Fmi 

1° 


max if ^ [ ^max5^max]i 

otherwise. 


We choose to present our results with the following parameters 

^max = 7 r/ 6 , a = 0.8 and L = 407tcj -1 . 


(23) 


(24) 


Because the evolution of different neutrino moments is decoupled in the linear regime, it is sufficient to include only 
the Oth and rath moments in studying the evolution of the rath moment in this regime. [The Oth moment is needed 
because it has large diagonal elements even in the linear regime. See Eq. @•1 

In Fig. [2] we show the numerical solutions to Eq. 0 in two calculations with all but the Oth and 1000th moments 
being zero. In both calculations, |e m ^|, the amplitudes of the off-diagonal elements of grow exponentially which 
is understood as flavor instabilities. As comparison we plot in Fig. [2] the exponential growth functions ~ exp(ft™ ax z) 
predicted by the flavor stability analysis, and they agree with the numerical results very well. As a further confirmation, 
we have compared the shapes of |Q m? $| and |Q m? $| obtained from flavor stability analysis (the dotted and dashed 
curves in Fig. [3| with those of |e m? #| and \e m ^\ in numerical calculations (not shown), and they also have good 
agreement. However, to achieve numerical convergence a large number of angle bins may be needed for the following 
reason. 

As pointed out in Ref. m there can exist many spurious flavor instabilities in the numerical implementation using 
the discrete (angle-bin) scheme. This can be seen from the discretized version of Eq. (10): 


i cos sin^j -f A cut] T (1 

- fj ,- cos(i?* - #j)](e mit?j . - ae m ' 0 j )g(‘&j), 
j 

icostiidzEmfo = [k m sin + X + LOT] + (1 - a)fi^ i }e rn ^ i 
^ ' [1 C0S($, 0j )] {^rn/Jj 


(25a) 


(25b) 


or 


idz^-m — An 


(26) 
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FIG. 3. The amplitudes of the unstable modes of the rath neutrino moments (in arbitrary scale) as functions of neutrino 
emission angle ft which have the largest exponential growth rates in the linear regime at given neutrino number densities 
(indicated by — ^/2Gfti u which is measured in the vacuum neutrino oscillation frequency u). The top and bottom panels 
are f or t he inverted and normal neutrino mass hierarchies, respectively. In these calculations we used the parameters listed in 
Eq. ( [24| , and we took the matter potential A = 0. 


where e m = (e m?1?1 , , e m?1?2 , e m ^ 2 ,..., e m ^ N , em,v N ) T is a 27V-dimensional vector, and A m is a 27V x 2 N real 

matrix. Matrix A m has 2 N eigenvalues Qm (i = 1 , 2 ,..., 2N) each of which corresponds to the collective oscillation 
frequency of a collective mode in the discrete scheme. Many of these collective modes can be unstable, i.e. with 
Km = Im(Q^) > 0. Only a few of the unstable modes correspond to the physical instabilities in the continuum limit 
(of the ft distribution), and the rest of them are “spurious” or the artifact of the numerical implementation. 

In Fig. [ 4 ] we plot the exponential growth rates Km of all the unstable collective modes both in the discrete scheme 
and in the continuum limit for the 0th and 5000th moments, respectively. This figure shows that spurious instabilities 
(in the discrete scheme) can dominate the physical instabilities (in the continuum limit) on small distance scales 
and/or large neutrino number densities (i.e. large \m\ and/or fi). In some extreme cases, e.g., the bottom middle 
panel of Fig. [4] where 77 = +1, m = 5000 and N = 100, none of the collective modes in the discrete scheme matches 
the ones in the continuum limit. This is likely due to the fact that Q m ,tf and Q m # become sharply peaked functions 
of ft at large \m\ and/or /r, which requires more angle bins to resolve (see Fig. |f|. Indeed, the comparison between 
the middle and right panels of Fig. [4] show that the spurious instabilities are more suppressed when more angle bins 
are employed. 


B. Flavor instabilities and matter effect 


We have solved the flavor instabilities of the multi-angle Line model using the angular distribution in Eq. (23) and 


the parameters listed in Eq. (24). The results for the neutrino gas in the absence of matter are shown in the upper 
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FIG. 4. The exponential growth rates of the unstable collective modes of the rath neutrino moment as functions of 
neutrino self-interaction strength /jl — \Z2Gfti u in the discrete angle-bin scheme with N angle bins (as labeled and shown as 
the dotted curves) and in the continuum limit of angular distribution (solid curves), respectively. Both k and /jl are measured 
in the vacuum neutrino oscillation frequency u. The top and bottom panels are for the inverted and normal neutrino mass 
hierarchies, respectively. In these calculations we used the parameters listed in Eq. (24), and we took the matter potential 
A = 0. 


panels of Fig. [5j From this figure one can see that, unlike the two-beam Line model El, the flavor instabilities in 
the multi-angle model depends on the neutrino mass hierarchy, and collective oscillations can begin at larger neutrino 
density in NH than in IH. One also sees that both /i™ ax and the maximum and minimum fi values where the 

mth modes are unstable, seem to increase linearly with \m\. In contrast, both /i™ ax and /i™ m increase linearly with 
VH in the two-beam model0 This implies that, for sufficiently large |m|, flavor instabilities can develop at even 
larger neutrino densities in the multi-angle model than in the two-beam model. 

Unlike in the two-beam model, the presence of matter can affect collective oscillations in the multi-angle model 
because the neutrinos propagate in different directions can travel through different distances between two lines that 
are parallel to the neutrino Line. In the lower panels of Fig. [5] we show the flavor instabilities in the multi-angle Line 
model with A = 20CL;. Similar to the situation in the spherical neutrino Bulb model for supernova mm, both /i™ ax 
and /i™ m of the homogeneous mode (i.e. with m = 0) shift to larger values in the presence of a large matter density in 
both NH and IH. However, /i™ of inhomogeneous modes actually shifts to smaller values for both NH and IH when 
\m\ is sufficiently large. 


IV. DISCUSSION 

We have used both the numerical method and the linear stability analysis to investigate collective neutrino os¬ 
cillations in the multi-angle Line model in the linear regime where the neutrino flavor transformation is still small. 
Although the Line model does not represent any real physical environment, the study of this toy model can provide 
insights into the important differences between the models of one spatial dimension (e.g. the neutrino Bulb model for 
supernova) and multi-dimension models. 


1 The definition of neutrino self-coupling strength /i in El, has taken into account of the geometric factor 1 — cos($ — $') and is equivalent 
to jltf in this paper. For the angular distribution in Eq. ( |23| jl o = (1 — sin7? max /'i9max)jw ~ 0.045 ji. 
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FIG. 5. Maximum exponential growth rate /t” ax (A, /x) (indicated by the color scale) of the neutrino collective oscillation modes 
in the multi-angle Line model as a function of moment index m and the neutrino self-coupling strength /x = \/2Gfti v . Both k 
and fi are measured in the vacuum neutrino oscillation frequency u. The left and right panels are for the inverted and normal 
neutrino mass hierarchies, respectively, and the top and bottom panels are for A = \f2GFTi e — 0 and 200a;, respectively. In 
these calculations we assume isotropic neutrino fluxes within angular range $ G [— 7t/6, 7t/6], and we used the parameters listed 
in Eq. d24}. 


An important goal of this work is to check if the inhomogeneous collective modes are suppressed in the multi-angle 
environment because of the high neutrino densities which is known to exist in the Bulb model nans]. Somewhat 
surprisingly, our work suggests that, in the absence of ordinary matter, inhomogeneous collective modes on small scales 
are not only not suppressed in the multi-angle environment, but can become unstable at larger neutrino densities 
than in the two-beam model. 

We also examined whether the presence of a large matter density can suppress collective oscillations in the two- 
dimensional Line model as in the one-dimensional Bulb model m- Our study shows that the presence of ambient 
matter does suppress inhomogeneous oscillation modes on large distance scales in the Line model as it occurs to the 
homogeneous modes in the Bulb model. However, it appears that the inhomogeneous modes on very small scales 
can occur at smaller neutrino number densities with ambient matter than without. In addition, the flavor unstable 
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region of the certain inhomogeneous modes can extend to the regime of lower neutrino densities than that for the 
homogeneous mode. 

We have shown that, as in the Bulb model, there exist spurious oscillations in the numerical implementation of the 
multi-angle Line model if the discrete angle-bin scheme is employed. The problem of spurious oscillations appear to 
be more severe at higher neutrino densities and on smaller distance scales. Although this problem can be mitigated 
by using more angle bins, it does add complications to the already challenging task of computing collective neutrino 
oscillations near astrophysical neutrino sources such as core-collapse supernovae and black-hole accretion discs. It is 
probably helpful to develop the multipole expansion method similar to that for the Bulb model m • 

Our work has focused on the neutrino flavor instabilities in the linear regime. However, not every flavor instability 
in the linear regime can result in significant neutrino flavor transformation. For example, in the realistic supernova 
environment, the neutrino density decreases as neutrinos travel away from the center of the supernova which results 
in the shift of the instability region. It is, therefore, possible that a collective oscillation mode does not grow all the 
way to the nonlinear regime during the finite distance interval where it is unstable. We have considered the mixing 
of two neutrino flavors only, which can be quite different from the neutrino flavor transformation of three flavors ED 
Ultimately, the phenomenon of collective neutrino oscillations has to be studied in realistic, multi-dimensional models 
for compact objects such as core-collapse supernovae and black-hole accretion discs before one can fully understand 
the impact of neutrino oscillations to these extreme environments. 
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